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In a recent paper Q] it was shown that a simple mean-field dynamo model with a spherically 
symmetric helical turbulence parameter a can exhibit a number of features which are typical for 
Earth's magnetic field reversals. In particular, the model produces asymmetric reversals (with a 
slow decay of the dipole of one polarity and a fast recreation of the dipole with opposite polarity), 
a positive correlation of field strength and interval length, and a bimodal field distribution. All 
these features are attributable to the magnetic field dynamics in the vicinity of an exceptional point 
of the spectrum of the non- self adjoint dynamo operator where two real eigenvalues coalesce and 
continue as a complex conjugated pair of eigenvalues. Usually, this exceptional point is associated 
with a nearby local maximum of the growth rate dependence on the magnetic Reynolds number. 
The negative slope of this curve between the local maximum and the exceptional point makes the 
system unstable and drives it to the exceptional point and beyond into the oscillatory branch where 
the sign change happens. A weakness of this reversal model is the apparent necessity to fine-tune the 
magnetic Reynolds number and/or the radial profile of a in order to adjust the operator spectrum 
in an appropriate way. In the present paper, it is shown that this fine-tuning is not necessary in 
the case of higher supercriticality of the dynamo. Numerical examples and physical arguments are 
compiled to show that, with increasing magnetic Reynolds number, there is strong tendency for the 
exceptional point and the associated local maximum to move close to the zero growth rate line where 
the indicated reversal scenario can be actualized. Although exemplified again by the spherically 
symmetric a 2 dynamo model, the main idea of this "self-tuning" mechanism of saturated dynamos 
into a reversal-prone state seems well transferable to other dynamos. As a consequence, reversing 
dynamos might be much more typical and may occur much more frequently in nature than what 
could be expected from a purely kinematic perspective. 

PACS numbers: 47.65.+a, 91.25.-r 



I. INTRODUCTION 

The magnetic field of the Earth is known to undergo 
irregular reversals of its dipole part. The reversal rate is 
variable in the course of time: it was nearly zero in the 
Kiaman and the Cretaceous superchrons and is approxi- 
mately 5 per Myr in the present Q ■ 

Much effort has been devoted to identify typical char- 
acteristics of the reversal process. In particular, it was 
claimed that reversals may have an asymmetric, saw- 
toothed shape, with the field of one polarity decaying 
slowly and recreating rapidly with opposite polarity, pos- 
sibly to quite high intensities 0, 0, HJ • 

Another hypothesis concerns a correlation between the 
polarity interval length and the magnetic field intensity 
|(| Q • Associated with this, a causal connection of the 
Cretaceous superplume and superchron period has been 
discussed d, |j| ■ It was Vogt who first suggested a cor- 
relation between volcanism and reversal rate 0|. The 
general idea behind this is that superplumes give rise to 
an increased heat transport from the core mantle bound- 
ary to the Earth surface with the result of an increased 
dynamo strength due to a higher temperature gradient 
driving the outer core flow 111. While this idea was 
soon generally accepted (with some counter-arguments 
regarding the involved time-scales [T^p. quite contrary 
implications for the reversal frequency were drawn from 
it. The first "school", advocating a negative correlation 
of interval length and energy supply to the dynamo, goes 
back to Loper and McCartney [13j. The second school, 



suspecting long intervals for a strong dynamo, was mo- 
tivated by various mean-field dynamo models, for which 
a transition from anharmonic oscillations to superchrons 
for increasing magnetic Reynolds number was observed 
0. 

A third, and still controversially discussed observation 
concerns the bimodal distribution of the Earth's virtual 
dipole moment (VDM) with two peaks at about 4 x 10 22 
Am 2 and at about twice that value 0, 0, . 

For decades, it has been a challenge for dynamo the- 
oreticians to explain reversals and their characteristics. 
It was considered a breakthrough when Glatzmaier and 
Roberts observed a reversal process in their fully coupled 
three-dimensional simulation of the geodynamo 18] (cf. 
also 0] for a recent overview). The strange thing with 
these simulations is that they reproduce many features 
of Earth magnetic fields, including reversals, quite well 
despite the fact that they are working in parameter re- 
gions far beyond those of the real Earth. This deficiency 
applies, in particular, to the Ekman and the magnetic 
Prandtl number. A way out of this dilemma may lie with 
a reliable sub-grid scale modeling 20] . In this respect one 
should also notice recent efforts to link direct numerical 
simulations and mean-field dynamo models 0, ■ 

This brings us back from expensive simulations to the 
complementary tradition of understanding reversals in 
terms of reduced dynamo models. A very simple ap- 
proach in this direction is the celebrated Rikitake dy- 
namo of two coupled disk dynamos [23L I24I] . 

Another model was studied by Hoyng and collabora- 
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tors [2^, [2|| • A mean-field dynamo model is reduced 
to an equation system for the amplitudes of the non- 
periodic axisymmetric dipole mode and for one periodic 
overtone under the influence of stochastic forcing. This 
simple model, which produces sudden reversals and a 
Poissonian distribution of the interval time, has also been 
employed to simulate the phenomenon of stochastic reso- 
nance |28j . Stochastic resonance was made responsible in 
a former paper 2 9| for an apparent 100 kyr periodicity in 
the interval length distribution l30j fnote however, that 
this periodicity is not settled yet |31|). An essential in- 
gredient of Hoyng's model to explain the correct reversal 
duration and the interval length consistently is the use 
of a large turbulent resistivity which is hardly justified. 
At least nothing of this has been seen in the recent liquid 
sodium dynamo experiments |3^| . 

A further approach to understand reversals relies on 
the transition between non-oscillatory and oscillatory 
eigenmodes of the dynamo operator [33l 133 . |3f| . Those 
transition poi nts, which have been found in many dy- 
namo models |36l l37| . are well known in operator the- 
ory as spectral branch points — " exceptional points" of 
branching type of non-selfadjoint operators |38| . Such 
branch points are characterized not only by coalescing 
eigenvalues but also by a coalescence of two or more (geo- 
metric) eigenvectors and the formation of a non-diagonal 
Jordan block structure with associated vectors (algebraic 
eigenvectors ] [3 9l Eel EH ] . This is in contrast to " diabol- 
ical points" [42| which are exceptional points of an acci- 
dential crossing of two or more spectral branches with an 
unchanged diagonal block structure of the operator and 
without coalescing eigenvectors [3|| E3 ■ 

In a recent paper 0], we have analyzed the magnetic 
field dynamics in the vicinity of an exceptional point in 
more detail. Although the used model, a mean-field dy- 
namo of the a 2 type with a supposed spherically symmet- 
ric helical turbulence paramter a, is certainly far beyond 
the reality of the Earth dynamo (owing, in particular, to 
the missing North-South-asymmetry of a) it exhibited all 
mentioned reversal features: asymmetry, a positive cor- 
relation of field strength and interval length, and bimodal 
field distribution. 

All those features together were attributed to the very 
peculiar magnetic field dynamics in the vicinity of an 
exceptional point. Usually, this exceptional point is as- 
sociated with a local maximum of the growth rate curve 
at a slightly lower magnetic Reynolds number. 

If this local maximum lies above zero, then there is a 
stable fixed point to the left of it. However, any prevailing 
noise can trigger the system to switch to the unstable 
fixed point at the right of the local maximum. From 
there the system is driven, in an self-accelerating way, 
to the exceptional point and beyond into the oscillatory 
branch where it undergoes the very polarity change and 
comes back to one of the fixed points. 

If the local maximum is below zero, the system un- 
dergoes an anharmonic oscillation (a limit cycle) with a 
pronounced asymmetry of the "reversal" . However, noise 



can lift the local maximum above zero making the sys- 
tem stay in the fixed point for a while before resuming 
the anharmonic oscillation. 

This reversal scenario seems not unrealistic since it 
exhibits at least three reversal characteristics with only 
demanding the existence of an exceptional point and a 
nearby local maximum of the growth rate. The bad news 
is that it seems to require an artificial fine-tuning of the 
intensity and/or the spatial distribution of the dynamo 
source, in order to position the exceptional point and its 
accompanying local maximum close to the zero line. By 
checking a variety of a profiles in the kinematic regime, 
we have indeed observed that the spectral structure as 
it is necessary for reversals to happen occurs only sel- 
dom. Hence, the criticism that our particular choice of 
a(r) has not enough geophysical background to explain 
reversal [2l| , seems well justified. 

Thus motivated, it is the main goal of the present paper 
to find better arguments why a dynamo operator should 
have a reversal-prone spectrum. A first hint on the solu- 
tion of this puzzle can be found in the papers by Bran- 
denburg et al. 01 an d by Meinel and Brandenburg [Sj . 
For a mean-field disc dynamo (which is rather a model 
for galactic than for planetary dynamos), it was shown 
|44j that a dynamo in the highly supercritical regime can 
exhibit a pronounced reversal behaviour, although this 
would not be expected from considering only the kine- 
matic profiles of a. 

In the present paper we will verify if this behaviour 
in the highly supercritical regime can be understood in 
terms of the " exceptional point model" . Our main out- 
come will be that there is a strong tendency of saturated 
dynamos to evolve into a reversal-prone state where the 
condition, that the exceptional point and the local maxi- 
mum are situated close to the zero line, is indeed fulfilled. 
Applied to the Earth, a possible conclusion could be that 
reversals are rather due to a self-tuning of the saturated 
dynamo operator than to an accidental fine-tuning of the 
kinematic dynamo operator. 



II. THE MODEL 

Instead of simulating the real Earth dynamo with a 
fully coupled three-dimensional solver, we employ a very 
simple mean-field dynamo model in order to work out 
clearly the basic mechanism of reversals. 

The considered a 2 dynamo with a spherically symmet- 
ric helical turbulence parameter a leads to a system of 
two coupled partial differential equations with only one 
spatial variable (the radius) . This model is simple enough 
to allow for long-time simulations providing reasonable 
reversal statistics, but at the same time it is still com- 
plex enough to catch the essence of /lydromagnetic dy- 
namos. In contrast to their technical counterparts, the 
saturation of hydromagnetic dynamos relies strongly on 
the deformability of the dynamo source which is, in our 
case, the variable radial dependence of a. Note that in 
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the quite different context of the Riga dynamo experi- 
ment a similar one-dimensional back-reaction model was 
used which reflects also the hydromagnetic character of 
this dynamo experiment in contrast to experiments with 
a more constraint geometric flexibility |45l l46j. 

The magnetic field evolution of a kinematic a 2 dynamo 
is governed by the induction equation 



^— = V x (aB) + AB 

OT /!o<7 



(1) 



with a denoting the helical turbulence parameter which 
may depend on the position r and the time r 47] . The 
dynamo acts, within a sphere of radius R, in a fluid with 
electrical conductivity a. The magnetic field has to be 
divergence-free, V ■ B = 0. In what follows, the length 
will be measured in units of R, the time in units of [i^oR 2 , 
and the parameter a in units of (//oci?) -1 . Note that for 
the Earth we get a typical time scale fioaR 2 ~ 200 kyr, 
which results in a free decay time of 20 kyr for the dipole 
field. 

It is convenient to decompose B into a poloidal and a 
toroidal field component according to B = — V x (r x 
V5) — r x VT. The defining scalars S and T are then 
expanded in spherical harmonics of degree I and order m 
with the expansion coefficients s; |fra (r,r) and ti lTn (r,T). 

For the present case with a(r) = a(r), the induction 
equation decouples for each I and m into the following 
pair of equations: 



dr 
dr 



- — (rsi) 
r dr 2 



1(1 + 1) 



si + a(r,T)ti, (2) 



l__d 

r dr 
1(1 + 1) 



^-(rti) - a(r,T)-^-(rs l ) 
dr dr 



[ti - a(r, t)si] 



(3) 



The boundary conditions are dsi/dr\ r= i + (I + l)s;(l) = 
ti(l) = 0. In the following we focus our attention on 
the dipole mode with 1 = 1, and will henceforth skip the 
corresponding subscript of s and t; s := s\ and t := t±. 
The absence of the order m in Eqs. (2,3) follows from 
the spherical symmetry of a. It implies, in particular, 
a complete degeneration of axial and equatorial dipole 
modes. It is clear that for any more realistic model (e.g. 
with inclusion of the North-South asymmetry of a) this 
degeneration would be lifted. 

Let us assume that the profile of a in the kinematic 
regime, akin(r), represents a supercritical dynamo. Af- 
ter self-excitation has occurred, magnetic field saturation 
is ensured by quenching the parameter a. We do this 
with the angularly averaged magnetic field energy which 
can be expressed in terms of s(r, t) and t(r, t) |4g. This 
averaging over the angles, which represents a severe sim- 
plification, has been introduced in order to remain within 
the framework of the spherically symmetric a 2 model. In 
reality, of course, any quenching would introduce terms 
breaking the spherically symmetry of a. 



In addition to the quenching effect, we assume the a- 
profile to be affected by " blobs" of noise which are consid- 
ered constant within a correlation time T corr . Physically, 
this noise is not unnatural: it could be understood as 
a consequence of changing boundary conditions for the 
core flow, but also as a shorthand for the omitted influ- 
ence of higher multipole modes on the dominant dipole 
mode. 

Putting all together, a(r) takes on the time dependent 
form 



a(r, t) = C 



akin (r) 



l + E mag (r,T)/E, 



u +S(r,r), (4) 

mag 



where E mag is the magnetic energy, averaged over the 
angles, 



E mag (r, t) 



2s 2 (r,r) , 1 fd(rs(r,r))\ 2 , a 



+ t 2 (r,r) .(5) 



In the numerical scheme, the noise term S(r, r) will be 
treated in form of a Taylor expansion, 



S(r, t) = & (t) +Ur)r 2 + & (r) r 3 + £ 4 (r) r 4 



(6) 



with the noise correlation given by (^i(T)^j(r + ti)) = 
D 2 (l - |ri|/r corr )6(l - |ti|/t co „.)(^ . 

In summary, our model is governed by four parameters, 
the magnetic Reynolds number C, the noise amplitude D, 
a mean magnetic energy i?„ ag in the saturated regime, 
and the noise correlation time T corr . 

The equation system (2)-(4) is time-stepped using an 
Adams-Bashforth method. For the following examples, 
the correlation time T corr has been set to 0.02, and E^g 
has been chosen to be 100. The details of these choices 
are not very relevant. Roughly speaking, a shorter corre- 
lation time T corr would require a stronger noise amplitude 
D in order to yield the same effect. 



III. ANHARMONIC OSCILLATIONS IN THE 
VICINITY OF EXCEPTIONAL POINTS 

In we had considered the following particular form 
of the kinematic a profile: 

a kln (r) = C(-21.5 + 426.4 r 2 - 806.7 r 3 + 392.3 r 4 ). (7) 

Actually, this strange-looking Taylor expansion was the 
result of an Evolution Strategy search for oscillatory 
spherically symmetric a 2 dynamos 49|. For this a(r) 
profile, the growth rate dependence on C is shown as 
curve Ki in Fig. lb. "K" stands for kinematic, and the 
subscript " 1" indicates the eigenfunction with the radial 
wavenumber 1. Correspondingly, K2 denotes the eigen- 
function with the radial wavenumber 2. Ki and K2 coa- 
lesce at the (leftmost) exceptional point E and continue 
as a pair of complex conjugate eigenvalues. At the upper 
exceptional point E', the curves Ki and K2 split off again 
into a pair of real eigenvalues. 
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The critical value for this a(r) profile, C = 1, lies on 
the oscillatory branch, enclosed between the two excep- 
tional points E and E'. For slightly supercritical values 
of C, the time evolution is a nearly harmonic oscillation 
that becomes more and more anharmonic and asymmet- 
ric with increasing C. 

At C = 1.2785 this oscillation acquires nearly a rect- 
angular form (Fig. la). For this extremely anharmonic 
oscillation we analyze, at the eight different instants 1...8 
shown in Fig. la, the corresponding instantaneous pro- 
files a(r,r) (Fig. 2a). Governed by Eqs. (4) and 
(5), these instantaneous profiles depend on the instan- 
taneous magnetic field variables s(r, t) and t(r, r), which 
are shown in Figs. 2b and 2c, respectively. For each 
of these instantaneous a profiles, the full squares on the 
vertical dashed line in Fig. lb show the resulting instan- 
taneous growth rate. In order to identify the position 
of these points with respect to the exceptional point, we 
put them in the context of a whole growth rate curve, by 
artificially replacing the value C = 1.2785 of the actual 
quenched a profile by values between 0.7 and 1.4. This 
way we obtain the thin lines 1...8 in Fig. lb which makes 
it possible to interprete a reversal in terms of the consec- 
utive deformation of the original kinematic growth rate 
curve. 

Let us begin with the instant 1. The magnetic field 
(Fig. 2b, 2c) is high, therefore the quenching of a is 
also quite strong (Fig. 2a). The resulting instantaneous 
growth rate (Fig. lb) sits nearby the local maximum 
which is slightly below zero. Hence, the field starts to 
decay slowly. 

At the instant 2, the field has already weakened, the 
quenching of a is less pronounced as before, and the point 
of the instantaneous growth rate has moved in between 
the local maximum and the exceptional point. At the 
instants 3 and 4, characterized by even weaker fields and 
less quenching of a, the growth rate point has reached 
the oscillatory branch. It follows a short intermezzo in 
the higher non-oscillatory branch (instant 5) where the 
a profile has nearly taken on the unquenched shape. Af- 
ter this, the system returns via the oscillatory branch 
(instant 6) to the reversed state at instants 7 and 8. It 
should be noted that this reversal scenario does not need 
any dramatic change or even a total sign reversal of a, 
as it was proposed e.g. in [5fj| . 

Readers familiar with the van der Pol oscillator [Hl| 
may notice that the signal form of Fig. la resembles 
strongly a "relaxation oscillation" (we thank Clement 
Narteau for drawing our attention to this point). In- 
deed, a closer inspection of the dynamics of the van der 
Pol oscillator shows a similar behaviour of the leading 
instantaneous eigenvalues. 

At C — 1.2789, the local maximum of the growth rate 
crosses the zero line, and instead of the considered limit 
cycle in form of an anharmonic oscillation we get two 
fixed points (the crossing points of the growth rate curve 
with zero). The fixed point at the left of the local maxi- 
mum is a stable one. So it will need some noise to jump 
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FIG. 1: (a) Part of an anharmonic oscillation showing the 
asymmetry of the reversals for the case C = 1.2785. The 
points 1...8 indicate the instants which are analyzed in the 
text, (b) Growth rates for the kinematic profile correspond- 
ing to Eq. (7) and for the quenched a profiles at the eight 
instants 1...8. The circle marks the critical point C = 1 for 
the kinematic a profile. "E" indicates, for each of the con- 
sidered a profiles, the lower exceptional point where the two 
first eigenmodes with radial wave numbers 1 and 2 coalesce. 
"E' " marks the second exceptional point, beyond which the 
two eigenvalues split and continue as real ones again. 



over the local maximum. In this regime, the interval 
length is not anymore governed by the intrinsic frequency 
of the anharmonic oscillation, but by the intensity of the 
noise which triggers transition from the stable fixed point 
to the unstable fixed point. 



IV. SATURATION INTO REVERSAL-PRONE 
STATES 

After the short review of Q given in the last section, we 
are left with the impression that the indicated reversal 
scenario depends heavily on an artificial fine-tuning of 
the shape and the intensity of the ct{r) profile. 

In this section we will try to find better arguments 
for dynamos to be in a reversal-prone state. For this 
purpose, we change our focus from slightly supercritical 
dynamos to highly supercritical dynamos. 

One comment is due in advance. The very particular 
a profile in Eq. (7) was the outcome of an Evolutionary 
Strategy search for an oscillatory and dominant dipole 
mode. In |49| it was shown that this double demand 
constrains the variety of possible a(r) profiles to a rather 
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FIG. 2: Instantaneous a profiles and denning scalars of the 
magnetic field at the eight instants indicated in Fig la. (a) 
a(r). (b) s(r). (c) t(r). 



thin corridor. It is much easier to find profiles a(r) with 
an oscillatory I = 1 mode for which some higher multipole 
modes with I = 2,3,.. are dominant. In order not to 
overcomplicate the problem, in the following discussion 
we lift the demand for a dominant dipole mode by simply 
omitting any consideration of higher dipole modes. 

Let us start with the classical profile a.kin{ r ) = C 
which is known to possess only real eigenvalues for all C 
|47j . For the sake of concreteness, we consider the criti- 
cal value C = 4.49 and the supercritical values C = 10, 
20, and 50. We solve the induction equation system (2) 
and (3) coupled to the quenching equations (4) and (5). 
The quenched a profiles, a(r) — C/(l + E mag (r) / ) , 
which are shown in Fig. 3b, are then multiplied by a 
scaling parameter C*, and the resulting instantaneous 
growth rate curves in dependence on C* are shown in 
Fig. 3a. As in the former section, this procedure is in- 
tended to identify the position of the actual growth rate 
point with respect to the exceptional point. 

Fig. 3b shows that the quenching of a is not homoge- 
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FIG. 3: (a) Growth rates for the profiles a(r) — C* ■ 
C / {I + Ema 3 {r) / E° mag ) with C =10, 20, 50. The asymptotic 
quenched state is in all cases non-oscillatory, (b) Kinematic 
and quenched profiles a(r) for the four considered values of 
C. 



nous along the radius. The saturation mechanism modi- 
fies the original constant a in such a way that there is a 
stronger suppression for smaller radii (evidently, because 
the magnetic field is strongest in the central part of the 
dynamo) . This modification of the shape of a has a re- 
markable consequence for the spectrum. In Fig. 3a we 
see that the growth rate curves acquire an exceptional 
point which is moving toward the zero growth rate line 
with increasing C . In this particular example, the excep- 
tional point does not drop below zero. Nevertheless, we 
will find later (Fig. 9a) that even in this case the noise 
can trigger reversals. 

The next example, akin = 5/3 C r 2 , is illustrated in 
Fig. 4 (the factor 5 /3 results from normalizing the radial 
average of a(r) to the value for constant a). The general 
tendency is the same as for the previous example akin — 
C, but now the exceptional point drops clearly below the 
zero line, leaving the maximum of this curve only slightly 
above zero. For C — 50, say, the dynamo will "sit" on 
the stable fixed point at the left of the local maximum. 
At this point the field is rather strong. However, any 
noise can trigger a move to the unstable fixed point at 
the right of the local maximum from where the reversal 
process can start. 

The next example is akin = 6/3 C r 3 . Here the local 
maximum drops below the zero line for C = 20 and 50 
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FIG. 5: Same as Fig. 3, but for a(r) = 6/3 C*Cr 3 /(l + 
E mag (r) / E^ag) for which the critical value is C = 7.84. The 
labels "20 (m)" and "50 (m)" refer to the maximally quenched 
a profiles during the anharmonic oscillation. 



FIG. 6: Same as Fig. 5, but for a(r) = 5.75/3 C* ■ C*(l - 
6r 2 + 5r 3 )/(l + E mag {r) / E^g) for which the critical value is 
C — 6.78, leading to an oscillatory mode. The label "7.2 (m)" 
indicates again the maximally quenched a. Note the move of 
the exceptional point well above the zero line and back to it. 



(Fig. 5). That means there is no stable fixed point any- 
more, the system runs into a limit cycle in form of the 
previously analyzed anharmonic oscillation. Therefore, 
a is changing its shape during this process. The curves 
"20 (m)" and "50 (m)" refer to those a profiles which are 
maximally quenched during the oscillation. 

The last example, afcm(r) = 5.75/3 C (1— 6r 2 + 5r 3 ), is 
quite similar to the example from but it will provide 
a surprise (Fig. 6). Starting from the kinematic a, for 
which the exceptional point is well below the zero line, 
it rises rapidly above zero to a maximum value, but for 
even higher C it moves back in direction of the zero line. 
Evidently, there are tow different mechanisms at work 
here. 

All four considered examples exhibit a tendency for su- 
percritical dynamos to saturate in a state for which the 
exceptional point and its associated local maximum lie 
close to the zero growth rate line. Interestingly enough, 
this happens independently on whether the exceptional 
point in the kinematic case was above the zero line (in- 
cluding the limiting case that there was no exceptional 
case at all) or below it. 

What is the physical rationale behind this phe- 
nomenon? Back-reaction for hydromagnetic dynamos is 
quite generally an actualization of Lenz's rule stating 
that the excited magnetic field acts against the source of 
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its own generation. For the considered model this means 
that the a profile is deformed in such a way that the 
growth rate is decreasing. In the first three examples 
this makes the growth rate curve of the leading eigen- 
mode come closer to the growth rate curve for the next 
but leading eigenmode (Fig. 7a). As a results, the back- 
reaction has a tendency to move the exceptional point 
downward in direction of the zero line (Fig. 7a). 

While this argument is rather intuitive for the first 
three examples, it does not apply to the fourth example 
in which an exceptional point is already existent below 
the zero line. Again, back-reaction is expected to lower 
the growth rate, but now the most efficient way to do this 
is quite different. Figure 7b may help to illustrate the es- 
sential point. The eigenvalues with radial wave numbers 
n = 1 and n = 2 have already coalesced and formed an 
oscillatory eigenmode below the zero line, but soon after 
they split again into a pair of real eigenvalues. The very 
steep increase of the leading eigenvalue in this upper real 
branch is quite typical and is the heritage of the level 
repulsion which would have occurred for a dynamo op- 
erator with a slightly modified a profile which is closer 
to the constant one. Therefore, the most efficient way 
of back-reaction in this case is to decrease the steepness 
of this leading eigenvalue branch. The decrease of the 
steepness goes hand in hand with a lengthening of the 
oscillatory branch. This saturation mechanism is so effi- 
cient that the local maximum and the exceptional point 
are even risen up towards the zero line and beyond. Only 
later, when the exceptional point is already well above 
the zero line, any further increase of C will promote the 
previous saturation mechanism, driving the exceptional 
point again back to the zero line (cf. Fig. 6). 

At present, work is in progress to support this rather 
intuitive and numerically derived picture by a paradig- 
matic analytical model. 

In some respect, this self-tuning mechanism is similar 
to the concept of self- organized criticality .52], although 
it seems too early to over-stress this resemblance. 



V. TIME SERIES AND REVERSAL DURATION 

The time evolution of the considered dynamos is, in the 
noise-free case, determined by the position of the local 
maximum of the growth rate curve relative to the zero 
line. If the local maximum is below zero, than we get 
a limit cycle in form of an anharmonic oscillation, the 
mechanism of which has been described in detail above 
and in . If the local maximum is situated well above 
the zero line, then the system runs into a stable fixed 
point. In case that the maximum is only slightly above 
zero, than a hysteretic behaviour may occur, in which 
the choice of the stable fixed point or the limit cycle is 
controlled by the strength of the pre-existing magnetic 
field. 

A quantity which is of particular geophysical relevance 
is the duration of a reversal. Here is not the place to 




FIG. 7: The self-tuning mechanism of reversing dynamos, (a) 
Case that there is no exceptional point or that it is above zero 
in the kinematic regime, (b) Case that the exceptional point 
is below zero in the kinematic regime. 



discuss the variety of different definitions, and all ques- 
tions concerning the dependence of the apparent reversal 
duration on the site latitude |53J. Since our model relies 
exclusively on the evolution of the dipole field, we cannot 
use any definitions based on directional changes. For that 
reason we had employed in Q a " working definition" of a 
reversal duration, based on the period during which the 
modulus of the magnetic field is smaller than a certain 
percentage of the amplitude of the (anharmonic) oscil- 
lation. For a value of 25 percent (which was motivated 
by demanding the dipole field to decay until this value 
in order that the non-dipole field can become dominant) 
we obtained a reversal duration of 0.15 Tdiff- For an as- 
sumed diffusion time of 200 kyr, this would amount to 
30 kyr. 

Now we would like to know how this reversal duration 
changes with increasing magnetic Reynolds number, in 
particular for highly supercritical dynamos. In Fig. 8 we 
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FIG. 8: Details of the reversal process for akin = 6/3 C r 3 
with the particular values C = 20, 50, 100, and 200. The 
time scale is the diffusion time, which is approximately 200 
kyr for the Earth. Evidently, the reversal duration decreases 
significantly with increasing C and reaches values of ~ 10 kyr 
for large C. 



plot the details of the reversal process for the third exam- 
ple of the previous section, aki n {r) = 6/3 C r 3 . For C we 
have chosen the four values 20, 50, 100, and 200. What- 
ever the exact definition of a reversal might be, the ten- 
dency is clear: with increasing C, the reversal duration 
decreases significantly. This has to do with the fact that 
during the reversal the magnetic field gets weak and the 
a profile comes close to its unquenched, kinematic shape 
for which extremely high instantaneous growth rates and 
frequencies (only in the oscillatory branch) occur. These 
are responsible for a very fast reversal process. If we 
would define the reversal duration as the period during 
which the field is between ± 50 percent of the oscillation 
amplitude, we would get for C = 200 a time span of 0.03 
T diff, corresponding to 6 kyr. The corresponding time 
for C = 20 is 0.09 r^//, i.e. 18 kyr. Evidently, we can 
get a realistic time scale without taking resort to turbu- 
lent resistivity, as it was done in other reversal models 

m 

So far, we have considered the noise-free case. As al- 
ready remarked, the role of noise depends on the position 
of the local maximum of the growth rate curve relative 
to the zero line. If it is above zero, the dynamo is usu- 
ally at the fixed point characterized by a strong magnetic 
field, then the noise can trigger transitions to the unsta- 
ble fixed point at the right of the local maximum, from 
where a reversal process can start. However, although 
being unstable, this fixed point can hold the system for 
a while (since the growth rate is zero there) making a 
second peak in the field strength histogram possible |rj. 
If the local maximum is below zero, making the dynamo 
undergo anharmonic oscillations with a weak field ampli- 
tude, then the noise can make it jump time by time to 
the strong field state. 

Anyway, the noise will soften the differences between 
the two regimes with the local maximum below or above 
zero. In either case, there will be noise-strength depen- 
dent exchange between the fixed point state with a strong 
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FIG. 9: Typical time series for the inclusion of noise, (a) 
a kln (r) = C. (b) a km {r) = 5/3 C r 2 (c) a ki „(r) = 6/3 C r 3 
(d) atkin(r) = 5.75/3 C (1 — 6r 2 + 5r 4 ) In all cases, we have 
chosen C = 50. In the case (a), we have set D = 7, in all 
other cases (b-d) D = 4. 



field and the limit cycle with a weak field. 

In Figs. 9 and 10 this is illustrated by some typical 
time series for the four dynamo examples considered in 
the previous section. In all cases we have chosen C = 50. 
The noise intensity, D has been set to 4, apart from the 
first example with akin — C for which a value of D = 7 
has been chosen to provoke any reversal at all. 

Apart from details concerning the mean reversal rate 
and the transient excursions to rather high values, all the 
signals share the property that there are long time inter- 
vals without reversals and rather fast reversal processes. 

A last remark concerns the long-standing problem of 
positive or negative correlation of energy supply and in- 
terval length. Our model is not able to solve this prob- 
lem, but it can illustrate its intricacy. The first point 
is that any additional energy supply will increase both 
the magnetic Reynolds number C and the noise level D. 
Even if considered separately from D, the influence of C 
depends on whether it will lift or lower the local maxi- 
mum with respect to the zero line. We have seen that 
both can happen, but for high supercriticality there is a 
general tendency of lowering the local maximum. This 
would imply a negative correlation of energy supply and 
interval length. If considering the influence of D sepa- 
rately from C, we get also a complicated behaviour with 
a tendency toward a negative correlation of D and the 
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FIG. 10: Same as Fig. 9, but zoomed into time periods were 
one reversal occurs. 

interval length, but with a positive correlation for the 
case that the local maximum is below zero and the noise 
is weak (cf. Fig. 7 of 0. Taken all together we get a 
very complex picture with no definite answer. 



reversals can a) be asymmetric, b) yield a positive cor- 
relation of dynamo strength and interval length, and c) 
show a pronounced bimodal field distribution. All three 
features have been recently discussed as being typical for 
reversals of the Earth's magnetic field. 

The apparent weakness of this model, the necessity of 
fine-tuning the spatial structure of the dynamo source 
in order to bring the exceptional point and its associ- 
ated local growth rate maximum close to the zero line, 
has been overcome in the present paper. We have shown 
that the back-reaction of the magnetic field has a strong 
tendency to drive the dynamo into a state where the indi- 
cated spectral conditions for reversals are indeed fulfilled. 

This shows that dynamos with a high, but not neces- 
sarily extreme, supercritical magnetic Reynolds number 
are very prone to reversals. Hence, the proposed reversal 
scenario which might look contrived and hardly realistic 
from a purely kinematic point of view, becomes rather 
natural when seen from the side of the saturation pro- 
cess. The artificial fine-tuning for the regime of slightly 
supercritical dynamos is replaced by a self-tuning satu- 
ration into reversal prone states when it comes to highly 
supercritical dynamos. 

The question remains if the Earth's dynamo is indeed 
highly supercritical. A first estimate of the magnetic 
Reynolds number, based on a length scale of 2000 km, 
a velocity scale of 5 mm/s, and a magnetic diffusivity 
of A := (nay 1 ~ 2 m 2 /s |H3, provides Rm ~ 500. At 
least, and apart from all uncertainties, this number does 
not exclude a highly supercritical state. Hence, whatever 
the concrete flow field in the Earth might be, it is not a 
surprise that the resulting dynamo is prone to reversals. 



VI. CONCLUSIONS 
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